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Measuring interface tensions in 4d SU(N) lattice gauge theories 
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We propose a new algorithm to compute the order-order interface tension in SU(N) lattice gauge theories. The 
algorithm is trivially generalizable to a variety of models, e.g., spin models. In the case N = 3, via the perfect 
wetting hypothesis, we can estimate the order-disorder interface tension. In the case N = 4, we study the ratio 
of dual k— tensions and find that it satisfies Casimir scaling down to T — 1.2 T c . 



1. Introduction 

The phenomenon of phase-coexistence is the 
typical signature of a first-order transition. The 
free energy of a system at the coexistence point is 
the sum of two contributions: the bulk free energy 
fbuikj scaling like the volume V, and the free en- 
ergy F-£ of the interface S separating the two bulk 
phases. scales like an area. The interface ten- 
sion a is defined as a = lims^co "^fsp (t ne reduced 
tension is defined as a R = §:)■ In the context of 
SU(N) pure gauge theories (JV is the number of 
colors) phase-coexistence occurs in two different 
situations: one at T > T c between ordered phases 
pointing in different directions in color space, the 
other at T c , between confined ('disordered') and 
deconfined ('ordered') phases. The correspond- 
ing interface tensions are indicated a QO and a d 
respectively. 

This paper presents a new algorithm to mea- 
sure the order-order interface tension. Neverthe- 
less, we can also give an estimate of the order- 
disorder tension via the so-called perfect wet- 
ting hypothesis pQ, which states that between 
two ordered phases a layer of disordered phase 
can be generated at no free energy cost. There- 
fore the free energy of one order-order interface 
is related to that of two order-disorder inter- 
faces, namely F ao — 2F 0( {. In general, anyway, 
Coo = wer d, w < 2, because otherwise order- 
order interfaces would be unstable. Therefore the 
choice w — 2 gives a lower bound for a 0( i- 



Figure 1. Cartoon of the interface, in an interme- 
diate situation in which only half of the plaquettes 
has been flipped. The solid central plaquette is 
affected by the least finite size effects. 



Our system is a Ad SU(N) pure gauge theory at 
finite temperature with spatial periodic boundary 
conditions (b.c). a indicates the lattice spacing. 
The volume is V = L 3 ■ L t ; the Wilson action is 
used. We will indicate the elements of the center 



of SU(N) by Ci, i 
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over-relaxation are applied to SU(2) subgroups 
of the SU(N) matrices, according to the strategy 
used in [5]. 



2. The method 

The algorithm we propose is an improvement of 
the so-called snake algorithm [3]. The 'snake' idea 
is to add 'by hand' a 2d interface in the system 
by progressively flipping the coupling of a set of 
plaquettes dual to a surface A, according to the 
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Figure 2. (Minus log of the) Ratios in Eq.QJ 
as a function of r, after removal of the Liischer 
correction. A broad plateau develops from r ~ £ 
(£ = -7= ~ 3.5 is the correlation length) to r — * L. 



identity 

Z(A) Z(A) 



Z(A-l) Z(l) 



Z(0) Z(A-l)Z(A-2) Z(0) 



(1) 



where Z(fc), fc <E [0, . . .A] indicates the partition 
function of a system in which only k plaquettes 
are flipped. The free energy of the interface is 
F(A) — cr^ Q A = — log^^j-. A direct measure- 



ment of Z z ^j is not possible due to a serious over- 
lap problem, which is alleviated by the factoriza- 
tion Eq . (|Tj) . The price to pay is a very large num- 
ber (L 2 ) of independent simulations. Progress 
in this direction is made comparing the following 
equations 



-oA 



Z(A) 
Z(0) 



n 



z(k) 

Z(k - 1) 



(2) 



the first just being a decomposition of the inter- 
face into the constituent plaquettes, the second 
indicating that each ratio, up to finite size (F.S.) 
corrections, contains all the information we need, 
namely 



Z(k-l) 



-(k— dependent F.S. corrections) (3) 



so that a single simulation suffices. The gain in 
efficiency is 0{L A ): a factor L 2 in the number 
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Figure 3. Finite size effects on the interface ten- 
sion (L t — 2, f3 — 5.15). After removal of the 
Liischer correction, a is systematically underesti- 
mated, unless lattices larger than L > 1 / \fa are 
used (vertical line). The curve is meant to guide 
the eye. 



of simulations, another factor L? because the L 2 
variances of the ratios add up in Eq.JQl. The 
observable is measured in the following way: 



Z{k) 
Z(k-l) 



( e /3±Tr(Cin fc )) fc 
/ e /3iTr(n fc )\ 



(4) 



where 11^ indicates the k— th flipped plaquette, 
and the average (•)& refers to the ensemble in 
which the first (k — 1) plaquettes are flipped, the 
A;— th plaquette has coupling zero, and all the oth- 
ers are unchanged. Further variance reduction 
methods are described in p{]. 

The leading finite-size corrections in Eq.© 
come from Gaussian fluctuations of the interface 
Our r ■ L interface (Fig.QJ is periodic in one 
direction (L), and pinned in the other (r). The 
corresponding correction, given in terms of the 
Dcdckind 77-function, reduces for r « L to the 
Luscher-like a e s(r) = -a~ 2 log Z(k)/Z(k — 1) w 
]S]. In Fig. |2] we show our measurements, 
after removal of this known correction, as a func- 
tion of r. A broad plateau develops, from small 
values r ~ £, where £ = l/v^ is the correlation 
length, to large ones r — > L, showing that addi- 
tional corrections are very small. At very large 
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Figure 4. Summary of results on the order- 
disorder interface tension for SU(3). To avoid 
an underestimate of o because of F.S. effects, we 
use a lattice of size 32 3 • 2 (a 2 T c 2 = 1/4), 48 3 • 3 
(a 2 T c 2 = 1/9), 64 3 ■ 4 (a 2 T 2 = 1/16). The latter 
simulation is in progress. 



distances, a systematic drop is visible, because it 
becomes more favorable to produce a full, trans- 
lationally invariant interface plus a partial one of 
width (L — r). Let us then fix r — (|f — 1); after 
removal of the Liischer corrections, Fig. [3] shows 
that the tension decreases considerably with the 
lattice size and reaches a plateau only when the 
empirical condition Lyfo > 7 is satisfied. 

3. Results 

In Fig. ^ we present preliminary results for 
the order-disorder interface tension in SU(3), to- 
gether with a compilation of the published data. 
While our Lt = 4 simulation needs more statis- 
tics, our Lt — 2 and 3 determinations of o are 
accurate, and much larger than previous mea- 
surements obtained with the histogram method 
EH 00 EH We assign this discrepancy to the 
smaller lattice sizes considered previously, which 
lead to a systematic underestimate of o as in 
Fig. El A discussion of the continuum limit is 
awaiting completion of the 64 3 • 4 simulations. 

In Tab n we present preliminary results for the 
case SU(4). In SU(N) with N > 3 we have 



T/T c 


L= 16 


L = 24 


2.3 


1.350(20) 


1.342(13) 


1.5 




1.300(18) 


1.2 


1.277(33) 


1.310(30) 



Table 1 

Ratio 02/01 of SU(4) interface tensions (L t = 5). 
The Casimir perturbative value is 4/3. For the 
data at T = 1.2 T c we find again that the F.S. 
effects are smaller than the statistical error only 



more order-order tensions (o&, k = 1, . . . , [N/2]). 
In the weak coupling regime one can show the 
Casimir relation 2i = k ^~^ 10 . This pertur- 
bative prediction seems accurate down to temper- 
ature T = 1.2 T c . Measurements closer to T c are 
in progress. 
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